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The general characteristics of stars with first-order phase transition (PT1) 

are reviewed at short. The model of two-incompressible-phase star with PT1 
is considered in some detail. 

I. INTRODUCTION 

First-order phase transition (PT1) is characterized by density jump from pi to P2 with 
Q = Pil P\ > 1) & t some pressure Pq. PT1 leads to density discontinuity inside the star while 
pressure is continious across boundary between phases. Also continious are gravity- and 
pressure-induced forces inside the star. Here we review some results of studying such stars 
get mainly (but not only) by author. 

II. THREE METHODS 

There are three main methods of star's equilibrium and stability analysis: 

• Static Criterion (Mass - Central Pressure dependence) 

• Dynamical Principle (pulsation frequency) 

• Variational (Energetic) Principle (variation of total energy). 

Stability loss/restoration (critical equilibrium states) according to these methods are 
denned as follows: 

• Mass extremum, dMj dP c = 

• Zero frequency of pulsation (of the lowest mode), uj 2 = 

• Equilibrium energy extremum, 5 2 E = 0. 
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In fig.l two first principles are shown qualitatively. Unstable equilibrium states are shown 
by broken line. 
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FIG. 1. Static (lower part) and dynamical (upper part) 
principles in the equilibrium and stability of stars 
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III. NEWTON THEORY OF GRAVITATION 



Here are some general results in NTG for stars with PT1: 

• At q > 1.5, stability loss occurs at P c = P for any EoS, while recover of stability for 
larger P c depends on EoS 

• At qminil) < 1 < 1-5, stability loss occurs only at P c > P ; 7 = 1 + l/n, 7 and n 
are adiabatic and polytropic indices. E.g., q min = 1.46, 1.33, 1.20, and 1.09 for n = 
1, 1.5, 2 and 2.5. 

The critical value of (relative) mass of new-phase core is larger for smaller g's and for 
smaller 7's (softer EoS's). 

E.g., critical value of new-phase core x crit = 8/9 (q — 3/2) for 7 = 1 (n — 1) 

• At slow rotation with angular velocity Q, in spherical approximation, q crit = 3/2 — 
VL 2 I 4 7r G pi, for stability loss at P c = P for any EoS. 

• For PT1, starting at finite radius ("neutral core"), 
the larger size of neutral core, the larger value of q crit . 
E~g-,C[crit — > 00 at x — > 3/4 at n = and x — > .6824 at n — 1. 

A. Classical example 

Inverse /3-decay reactions in dense degenerate matter of white dwarf stars lead to nuclei 
transformations (A, Z) — > (A, Z — 2) ("neutronization"), and to density jump with q = 
Z/(Z — 2) < 1.5. M — p c curves for equilibrium cold white dwarfs in classical paper by 
T.Hamada, E.Salpeter, Ap.J. 134 (1961) 669 are incorrect, as: 

a) central density p c is not continious variable, and b) mass maximum is not at point P c = P$, 
but at some P c > P . 
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B. General Relativity 



In GR, the critical value of energy density jump q = e 2 /£i is equal to 3/2 * (1 + Po/si) 
for stability loss at P c = P Q for any EoS. 

IV. TWO-INCOMPRESSIBLE-PHASE STAR 

Here are some important features of this model: 

• At q < 1.5 there is no unstable equilibrium states 

• At q > 1.5 recover of stability occurs (at point P c = P 2 ) for relative radius of new-phase 
core, x = r core /R sta r, defined by relation: 

f(q, x) = (q- 1) V + 4 (q - 1) x + 3 - 2q = (1) 

Star with x > ^2 — 1 is stable for arbitrary large q 

• Frequency squared of the small adiabatic radial pulsations of the lowest mode for star 
in slow rotation with angular velocity Q: 

= UJ% + A n (q,x); 

- 3( 9 - !)(!-*)• (2) 

At small cores (x — > 0), is negative - rotation reduces the stability of star with PT1. 
In general, Aq may be of both signs, e.g., at q > 2.11, rotation leads to decreasing of 
value of x cr u from Eq.(l). 
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A. Non-linear pulsations 



In the next approximaion, pulsations are non-harmonic. 
Writing down R = R eq + z, \z\ << R eq , R eq being radius of equilibrium model with the 
same mass, we get in next-to-zeroth approximation the following equation of motion: 

z + u%z + Cz 2 + Dz 2 = 0, (4) 

where u 2 is as in Eq. (2), while constants C and D are some functions of q and of equilibrium 
value of x. The solution of Eq. (4) with accuracy up to a 2 , ( a being an amplitude of z) is 
as follows: 

1 C 1 C 

z(t) = --(— + D)a 2 + a x cosw t + ^(— - D)a 2 + a 2 x cos2w £- (5) 

In this approximation the frequency does not differ from one in zero'th approximation, u Q , 
that is period T = T = 2ix/uq, while pulsations are non-sinusoidal. 

An amplitude of a star's expansion is larger than an amplitude of a star's compression 

and star spends more time with R > R eq than in state with R < R eq . A character of 

non-harmonicity - a slow large-amplitude "expansion" and rapid "contraction" with smaller 

amplitude - should be the same for all equilibrium stable models at A > 3/2. 

Damping of pulsations of star with PT1 depends largely on relation between velocity of 

motion v p and a sound velocity v s in the region of phase transition. 

In general, the larger v p /v s the larger damping effect due to PT. 

B. General Relativity 

Full analytical investigation of equlibrium and stability of two-imcompressible-phase 
model with PT1 is possible by static and variational methods. 

In first post-Newtonian approximation, (Pq/si 1) critical value of relative "radius" of 
core, at which stability recover occurs, 
is equal to: 
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XcriM) = X + A PN * (PoM), (6) 

with x = Xn defined in Eq. (1) and : 



9 - 7q + 27(q - l)x + (4g 2 - 33g + 27)x 2 + (g - 1)(9 - 4q)x 3 



±PN 



(7) 



2(g-l)(l + (g-l):r 3 ) 3 

At g -> 3/2, xpjv -> g - 3/2 - 3/2 P /ei, so that x P7V = at g = 3/2 (1 + P /£i), 
which coincides with result for am/ EoS and for any, not only for small values of Po/ei- In 
general, first post-Newtonian correction to x cr u may be of both signs, negative at q < 1.89 
and positive for larger q's. At q = 4 + y/8, x = ^2 - 1 + (59 - 41 v /2)/2 P /£i- 

In fact, dependence of both corrections, due to GR and rotation, on x and q is rather 
complicated, and Fig. 2 presents only part of (q — x)— plane with lines, on which A n (q, x) — 
and ApAr(g, x) = 0. Also shown is the curve f(q,x) = from Eq. (1) which marks a 
boundary between stable equilibrium states (at right) from unstable ones (at left). 




FIG. 2. Three important curves at (q — x)— plane: critical state of stability recover, see Eq. 1, 
and GR and Rotational corrections equal to zero, see Eq. (3) and Eq. (7) respectively. 
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